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摘 要 : 遥感 降水 产品 为 稀缺 资料 地 区 水 文 气象 灾 


害 机 理 与 预警 研究 提供 了 重要 的 数据 资料 ,但 不 


同 遥 感 降 水 产品 的 性 能 存在 较 大 区 域 异 质 性 。 利 用 遥感 降水 产品 开展 水 文 气象 相关 的 研究 和 应 
用 前 ,需要 对 其 性 能 进行 综合 评估 。 基 于 此 ,以 稀缺 资料 的 黄河 源 区 为 研究 区 ,利用 1983 一 2018 年 
的 观测 降水 数据 (CMA) 驱 动 并 率 定 ABCD 水 文 模型 ,并 利用 标准 化 径流 指数 (SRI) ,评估 3 套 典 型 
遥感 降水 产品 (PERSIANN-CDR .CHIRPS v2.0, MSWEP v2.0) 对 水 文 干旱 的 模拟 性 能 。 利 用 游程 理 
论 识别 水 文 于 旱 事件 ,阐明 遥感 降水 对 水 文王 旱 特征 的 捕 提 能力。 结果 表明 :(1)3 套 遥感 降水 产 
品 均 能 较 好 地 捕获 CMA 多 年 均值 的 时 空 分 布 格局 。CHIRPS 产品 的 水 文 模 拟 性 能 (纳什 效率 系数 
NSE=0.72) 高 于 其 他 2 套 产 品 ,(2) CMA 和 降水 产品 模拟 的 4 个 尺度 的 SRI(SRIl1、SRI3、SRI6 和 和 
SRI12) 均 呈 显 著 增 加 趋势 (P<0.01) ,表明 近 36a 源 区 河川 径流 增加 ,水 文 于 旱 趋 缓 , 但 降水 产品 均 


高 估 了 SRI, 表 明 对 黄河 源 区 降水 产品 的 偏差 校正 有 待 


展 。 基 本 统计 指标 方面 ,MSWEP 产品 计 


算 的 SRI 与 CMA 的 最 为 一 致 ,性 能 最 佳 ,但 在 年 尺度 (SRI12) 上 ,PERSIANN 表现 最 优 .(3) 3 套 产 品 
349 th Y SRI T SRB 的 干旱 历时 和 烈度 ,MSWEP 产 品 对 SRI6 的 模拟 性 能 最 优 ,PERSIANN 对 
SRI12 的 模拟 性 能 最 优 。 研 究 结 果 可 为 黄河 源 区 水 文 干旱 研究 的 降水 产品 数据 的 选择 提供 科学 决 


策 支 持 。 
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文章 编号 : 


在 气候 变化 和 人 类 活动 的 双重 影响 下 ,区域 
旋 至 全 球 的 平均 温度 和 极端 温度 均 呈 显著 增加 趋 
势 "。 温 度 的 急剧 升 高 加 快 了 水 循环 ,导致 干旱 事 
件 发 生 的 频率 、 历 时 和 强度 愈演愈烈 ,已 严重 威胁 
到 了 生态 文明 建设 和 经 济 高 质量 发 展 。 水 文 干旱 
作为 典型 干旱 类 型 之 一 ”, 是 气象 干旱 向 农业 干旱 
和 社会 经 济 干旱 传播 的 关键 过 程 ,研究 水 文王 旱 的 
演变 过 程 和 传播 机 制 有 利于 干旱 的 预警 监测 和 
防治 。 

为 了 量化 水 文 干 旱 程度 ,较为 常用 的 方法 就 是 
构建 水 文 干旱 指数 来 表征 ”。 在 众多 的 干旱 指标 研 
究 中 , 丁 品 等 “用 负 轮 长 指标 定量 研究 水 文 干 旱 , 分 
析 了 中 国 主要 河流 177 站 的 干旱 特性 ,但 该 方法 的 
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和 降水 构建 $ 指 数 , 并 评估 了 黄河 中 游 的 水 文 干 旱 
特征 ,但 是 对 $ 指 数 权 重 系数 的 确定 具有 较 大 难度 ; 
而 标准 化 径流 指数 (SRI)5 ,由 于 仅 需 要 径流 作为 输 
和 ,并且 计 算 方法 简便 和 计算 尺度 灵活 等 优点 , 现 
已 被 广泛 应 用 于 水 文 干旱 研究 之 中 。 为 表征 水 文 
干旱 的 演变 规律 和 驱动 机 制 ,和 常 采用 气象 数据 (如 
降水 . 洪 在 蒸 散 发 等 ) 驱 动 水 文 模 型 模拟 径流 过 程 ， 
进而 计算 SRI, 最 后 通过 数值 模拟 量化 不 同 因素 对 
于 SRI 的 影响 。 

降水 是 驱动 水 文 模型 的 核心 因素 ,其 精确 性 与 
否 会 给 干旱 事件 的 识别 带 来 很 大 不 确定 性 。 对 于 
降水 而 言 ,获取 方式 较 多 。 目 前 ,降水 数据 的 获取 
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主要 有 3 种 :地 面 观 测 .雷达 监测 和 遥感 卫星 反 演 ”。 
地 面 降水 观测 虽然 能 获得 降水 近似 的 真 值 ,但 由 于 
雨量 站 分 布 不 均 , 在 一 些 经 济 欠 发 达 地 区 和 地 形 复 
杂 区 ,雨量 计较 少 , 难 以 捕捉 降水 时 空 分 布 趾 ;天 气 
雷达 虽然 能 在 较 大 空间 范围 内 提供 较为 准确 的 降 
水 ,但 由 于 成 本 较 高 , 且 雷 达 容 易 受 到 多 种 误差 来 源 
影响 ,在 复杂 地 形 区 效果 不 佳 ”"; 相 比 之 下 遥感 降水 
产品 有 能 够 提供 精度 较 高 .时 间 序 列 连 续 的 降水 数 
据 的 优点 ,被 广泛 应 用 于 水 文 气象 的 研究 之 中 心 。 
目前 ,全球 遥感 降水 产品 众多 ,常见 的 如 PERSI- 
ANN-CDR, CHIRPS v2.0、MSWEP v2.0、CMORPH 
v1.0, GPCP-1DD v1.2, TMPA 3B42V7 等 i。 为 表征 
这 些 遥 感 降 水 产品 在 水 文 气象 领域 的 性 能 ,国内 外 
学 者 开展 了 评估 遥感 降水 产品 性 能 的 一 系列 研究 
工作 。 比 如 , 刘 洁 等 “采用 国家 基准 气象 站 观测 的 
降水 量 为 观测 数据 ,评估 了 TRMM 3B42RT TRMM 
3B42V7,CMORPH RAW 和 CMORPH CRT 4 套 遥 感 
降水 产品 在 塔里木 河流 域 的 性 能 ,发 现 TRMM 
3B42V7 最 适 于 模拟 径流 ;Bai 等 评估 了 CHIRPS 
V2.0, CMORPH V1.0, PERSIANN- CDR , TMPA 3B42V7 
All MSWEP V2.0 TE ji fi ibt E] FT X ed JY PERE , as 
AR aH] MSWEP 产 品 表现 最 优 ; 李 艳 忠 等 六 评估 了 
CHIRPS V2.0, CMORPH V1.0, TRMM 3B42, PERSI- 
ANN-CDR 和 MSWEPV2.0 产 品 在 渭河 流域 的 4 个 子 
流域 的 水 文 性 能 ,结果 表明 TRMM 的 水 文 模拟 性 能 
最 优 ; Guo 等 号 评估 了 PERSIANN-CDR , MSWEP 
V2.0 fll CHIRPS V2.0 在 西北 干旱 区 的 性 能 ,研究 发 
现 3 套 产品 在 估算 降水 和 表征 干旱 事件 方面 表现 
出 普遍 较 优 的 性 能 ; 王 蔡 等 "利用 TRMM 数据 和 
气象 站 点 数据 ,计算 标准 化 降水 指数 (SPI) ,分 析 了 
三 江 源 区 8 月 的 干旱 特征 。 综 上 可 知 , 较 多 研究 侧 
重 于 对 遥感 降水 产品 捕获 降水 的 时 空格 局 .径流 模 
拟 性 能 和 捕获 气象 干旱 性 能 的 分 析 ,而 对 遥感 降水 
捕获 水 文 干旱 性 能 的 研究 相对 较 少 。 此 外 ,不 同 的 
遥感 产品 在 不 同 地 区 存在 显著 的 区 域 差 异性 。 因 
此 ,开展 遥感 降水 产品 在 特定 的 区 域 捕获 水 文 干 旱 
的 性 能 ,有 助 于 拓展 遥感 降水 在 水 文 干旱 领域 的 
应 用 。 

黄河 源 区 是 黄河 的 发 源 地 ,面积 约 占 整 个 黄河 
流域 的 135%。 源 区 径流 量 约 占 黄河 径流 量 的 1/3 以 
上 ,车 有 ”黄河 水 塔 "之 称 ”"。 黄 河源 区 生态 环境 脆 
38 , 防 灾 减 灾 能 力 较 弱 后 , 且 对 气候 变化 和 人 类 活 


动 等 干扰 敏感 。 比 如 ,20 世 纪 90 年 代 , 由 于 气候 变 
化 和 过 度 放 牧 等 影响 ,以致 源 区 河川 径流 呈 显 著 下 
降 趋 势 ”” ,导致 了 水 文 干旱 的 发 生 。 稀 玻 的 降水 观 
测 站 点 数据 ,严重 制约 了 源 区 水 文 干旱 演变 过 程 和 
驱动 机 制 的 研究 ,而 遥感 降水 产品 为 开展 此 类 研究 
提供 了 可 能 。 面 对 诸多 的 遥感 降水 产品 后 , 哪 一 种 
产品 在 黄河 源 区 的 性 能 更 好 ? 哪 种 遥感 降水 产品 
更 能 准确 地 模拟 黄河 源 区 的 水 文 干旱 过 程 及 其 干 
旱 特征 ? 基于 此 ,本 研究 拟 选 取 3 套 典型 遥感 降水 
产品 (PERSIANN-CDR ,CHIRPS v2.0 .MSWEP v2.0), 
All FH ABCD 水文 模 型 ,探究 3 套 遥 感 降水 产品 在 黄 
河源 区 的 水 文王 旱 监测 性 能 ,为 黄河 源 区 水 文 干旱 
研究 的 降水 产品 数据 的 选择 提供 科学 决策 支持 。 


1 研究 区 概况 


BE Yi WRK (32°06' ~36° 12’ N, 95°48’ -103?24' E) 
位 于 青藏 高 原 东北 部 ,地 势 西高 东 低 , 源 区 流域 面 
FRIK 12x10* km (E 1). ZFFE 508 mm, 
多 年 平均 气温 -4.9 C ,平均 海拔 在 4000 m 以上， 
属于 高 原 亚 寒带 半 湿 润 气候 2。 源 区 内 植被 类 型 以 
高 寒 草包 为 主 ,区 域内 存在 季节 性 和 多 年 冻 土 区 。 


2 数据 与 方法 


2.1 数据 来 源 

准确 的 降水 数据 是 水 文 气象 研究 的 基础 ,涉及 
的 2462 个 站 点 的 逐日 气象 数据 来 自 国家 气象 信息 
中 心 (http://data.cma.cn/en) ,经 过 Anusplin 软件 插值 
形成 空间 分 辨 率 为 0.05°x0.05° 的 数据 ,时 间 范 围 为 


图 1 黄河 源 区 示意 图 
Fig. 1 Schematic diagram of the source region of 
the Yellow River 
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1980 一 2018 年。 由 于 Anuspline 方 法 考虑 到 高 程 和 
地 形 等 因素 ,相对 于 其 他 插值 结果 精度 更 高 下 。 唐 
JIZ Z 1980—2018 年 逐 月 径流 数据 来 自 中 国 水 利 
部 编 舞 的 水 文 年 鉴 。 本 研究 选取 了 3 套 典 型 遥感 降 
水 产品 ,包括 PERSIANN-CDR , CHIRPS v2.0, MSWEP 
v2.0( 以 下 简称 PERSIANN .MSWEP.,CHIRPS) 。 

PERSIANN 产品 由 Gridded 卫星 (GridSat-B1) 的 
红外 辐射 数据 作为 输入 数据 ,并 且 通 过 GPCP 降 水 产 
品 数据 进行 偏差 校正 。 该 产品 的 空间 覆盖 范围 为 
60?S~60*N2。CHIRPS 是 由 美国 地 质 调查 局 (USGS ) 
和 加 利 福 尼 亚 大 学 的 气候 灾害 小 组 共同 研发 的 一 种 
产品 , 它 融 合 了 地 面 雨量 站 数据 .卫星 数据 和 降水 气 
候 学 数据 等 ,空间 覆盖 范围 为 50°S~50°N。 

多 元 加 权 集 成 降水 MSWEP 是 一 套 全 球 降水 数 
据 集 ,其 独特 之 处 是 它 融 合 了 地 面 测量 、 卫 星 和 再 分 
析 数 据 。MSWEP 与 其 他 21 个 全 球 降水 产品 相 比 有 
着 更 加 卓越 的 性 能 ,并 且 在 缺 资料 或 无 资料 地 区 
具有 较 大 潜力 了。 为 了 评估 的 一 致 性 ,本 研究 空间 
4] BESSIS] Jy 0.25?x0.25? ,选取 的 研究 时 段 为 1983 一 
2018 年 。3 种 遥感 降水 产品 的 汇总 如 表 1 所 示 。 

2.2 研究 方法 

2.2.1 标准 化 径流 指数 (SRI) 标准 化 降水 指数 
(SPI) 是 评价 气象 干旱 的 常见 指标 。SPI 仅 基于 降 
水 序列 数据 计算 ,通常 选用 Camma 分 布 来 描述 降水 
量 的 变化 情况 ,能 够 将 偏 态 的 降水 分 布 转化 成 正 态 
分 布 进行 处 理 。SRI 是 评价 水 文 干旱 的 常见 指标 ， 
与 SPI 类 似 ,SRI 计 算 方法 只 需 将 降水 序列 数据 替换 
成 径流 序列 数据 , 石 朋 等 ”的 研究 表明 广义 极 值 分 
布 (GEV) 是 黄河 源 区 最 适合 的 流量 分 布 函 数 。SRI 
可 用 灵活 的 时 间 尺 度 来 量化 干旱 事件 ,考虑 了 1 个 
月 .3 个 月 .6 个 月 .12 个 月 4 个 典型 的 时 间 尺 度 的 
SRI(SRII SRI6 .SRI6 SRI12) 。 干 旱 等 级 划分 如 表 
2 所 示 。 

2.2.2 ABCD 水 文 模型 在 众多 水 文 模型 中 ,Thom- 
as 提 出 的 ABCD 水文 模型 有 参数 少 .模拟 精度 高 . 适 
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应 性 较 强 等 优点 ”。 模 型 的 输入 数据 为 降水 量 
(P) ,潜在 蒸 散 发 (PET) ,输出 数据 为 径流 深 (R) E 
壤 售 水 量 (S) 和 地 下 含水 量 (G)。 其 中 ,潜在 燕 散 发 
的 计算 采用 联合 国 粮农 组 织 开 发 的 Penman-Monte- 
it 公式, 计算 公式 如 下 : 

0.408(R, — G) + Y Dor udo, -e,) 


PELS A+y(1+0.34u,) a 
Wb W,-bY bW, 

n- 2a | 2a js a (2) 

R.=( - 9(W, - Y.) (3) 

R,-c(W, - Y.) (4) 

R, - dG, (5) 


式 中 :a 为 土壤 完全 饱和 前 形成 径流 的 概率 , 取 值 范 
围 为 (0, 11;5 为 不 饱和 含水 层 储 水 量 的 上 限 (mm)， 
取 值 范围 为 (0, 1000];c 为 土壤 含水 层 补给 地 下 水 
量 的 比例 , 取 值 范围 为 (0, 11;4 为 地 下 水 形成 出 流 
的 速度 , 取 值 范围 为 (0, 1j;R, 为 地 面 净 辐射 (MJ:m > 
d); Gg EROR EE (MJ -m7- d 2 ;y 为 干 湿 计 常数 
(kPa 5C) ;7 为 日 平均 温度 (%C) 5 us 23 2 m ig E Ab XU 
速 (m…s 1!);e 为 饱和 水 汽 压 (kPa) ;e, 为 实际 水 汽 压 
(kPa) ;A 为 水 汽 压 曲线 斜率 (kPa 9C) ;模型 的 2 个 
状态 变量 系 和 了 琴 分 别 为 可 能 薰 散发 量 (mm) 和 有 效 
7K HE (mm) ;RR. 为 直接 径流 (mm) ;RR 为 地 下 水 补给 量 
(mm) ;RR, 为 地 下 径流 (mm);G 为 当前 月 份 的 地 下 水 
储量 (mm )。 
2.3 统计 评价 指标 

为 定量 研究 3 套 遥 感 产品 在 黄河 源 区 的 水 文 干 
旱 性 能 ,选取 了 相关 系数 (CC ) 7T HEIR (RMSE) , 
Kling-Gupta 系数 (KGE) 纳什 效率 系数 (NSE)4 种 统 
计 评 价 指标 ,各 指标 的 计算 公式 如 下 : 


n 


»(6,- G)(S,-$ 
CC= HE“ NS (6) 


n 


[&-6r36.-3) 


n 
i= i=l 


(7) 


表 1 遥感 降水 产品 基本 信息 


Tab.1 Basic information of remote sensing precipitation products 


数据 集 年 份 空间 跨度 数据 下 载 来 源 
PERSIANN-CDR 1983—2018 60°S~60°N http://chrsdata.eng.uci.edu/ 
CHIRPS v2.0 1983—2018 50°S~50°N ftp://ftp.chg.ucsb.edu/pub/org/chg/products/CHIRPS-2.0 
MSWEP v2.0 1979—2018 全 球 http://gloh2o.org/ 
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表 2 标准 化 径流 指数 (SRI) 等 级 划分 标准 


Tab.2 Classification criteria of SRI 


—0.5<SRI 无 早 
-1.0<SRI<-0.5 轻 度 干旱 
—1.5<SRI<-1.0 中 度 干 旱 
-2.0<SRI<=1.5 严重 干旱 

SRI<-2.0 极端 干旱 


KGE-1- (1-CC)'+(1-a)' «(1-8Y, 
(8) 


式 中 : C. G 分 别 为 观测 降水 的 估算 降水 及 其 均值 
(mm); S, , S 分 别 为 遥感 产品 的 估算 降水 及 其 均值 
(mm) ;n 为 样本 容量 ;ow 分 别 为 观测 降水 的 均值 
和 标准 差 (mm) ;ow 分 别 为 遥感 降水 估计 值 的 均 
值 和 标准 差 (mm);Q.、8; 分 别 为 观测 径流 量 和 模拟 
径流 量 (mm); Q, 为 观测 径流 量 的 均值 (mm )。 


(a) CMA 


2.4 游程 理论 
干旱 事件 通过 游程 理论 确定 ,根据 McKee Ee 
的 研究 , 当 SRI 为 负 值 ,峰值 小 于 -1 且 干 旱 历时 持续 
2 个 月 以 上 时 ,就 被 认定 为 一 次 干旱 事件 ,记录 干旱 
事件 数 日 (DEN)。 干 量 历 时 (DD) 定 义 为 干旱 事件 
开始 和 结束 之 间 的 月 数 , 干 旱 烈度 (DS) 为 干旱 事件 
期 间 累 积 SRI 的 绝对 值 。 此 外 ,本 研究 选用 平均 干 
旱 历 时 (MDD) 和 平均 干旱 烈度 (MDS) 作 为 干旱 特 
征 指标 ,计算 公式 如 下 : 
wpp- LYpp, (10) 


1 N DD 
MDS = 4- È DS, DS= 2|SRI,| (11) 


式 中 :NN 为 干旱 事件 数量 ;DD 为 i 月 的 干旱 历时 ;DS 
为 干旱 事件 j 的 干旱 烈度 ;SRI 为 i 月 的 标准 化 径流 
指数 。 


3 结果 与 分 析 


3.4 遥感 降水 产品 时 空 分 布 
图 2 和 图 3 分 别 展 示 了 观测 降水 数据 (CCMA) 和 


(b) PERSIANN 


(c) MSWEP 


22 28 34 40 46 


ez 
ie 


Ps 


降水 量 /mm 
58 64 70 76 82 88 


图 2 黄河 源 观 测 降 水 数据 (CMA) 与 3 套 遥 感 降水 产品 的 月 均 降 水 量 空间 分 布 


Fig. 2 Spatial distributions of monthly average precipitation of CMA and three remote sensing precipitation 


products at the source of the Yellow River 


202308.00019v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


3 套 遥 感 降 水 产品 的 月 平均 空间 分 布 和 散 点 图 。 由 
图 可 知 ,3 套 遥感 降水 产品 均 能 够 捕捉 到 降水 量 自 
东南 向 西北 递减 的 空间 分 布 格局 。 但 3 套 遥 感 降 水 
产品 均 低 估 了 西北 部 的 降水 量 ,PERSIANN 产品 高 
估 了 东南 部 的 降水 量 。 此 外 ,3 套 遥 感 产 品 均 能 较 好 
地 捕获 到 CMA 的 波动 (CC>0.95) ,整体 CHIRPS 产 品 
的 性 能 表现 最 优 (KGE=0.95,RMSE=8.9 mm) ,其 次 
为 PERSIANN 产品 (KGE=0.94, RMSE=9.8 mm) , 较 
差 的 是 MSWEP 产 品 (KGE=0.95,RMSE=14.1 mm). 
3.2 ABCD 水 文 模型 性 能 

ABCD 水文 模型 率 定 期 和 验证 期 的 模拟 性 能 如 
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图 4a 所 示 ,模型 率 定期 (1980 一 2000 年 ) 的 NSE 达到 
0.81 ,验证 期 (2001 一 2018 年 ) 的 NSE 达 到 0.70 ,说 明 
模型 拟 合 较 优 , 表 明 ABCD 水 文 模型 在 黄河 源 区 有 
较 好 的 适应 性 ,可 用 于 不 同和 遥感 降水 产品 的 水 文 性 
能 验证 。 

基于 率 定好 的 水 文 模型 参数 ,将 3 套 遥 感 降水 
产品 驱动 ABCD 水 文 模型 ,其 模拟 性 能 如 图 4b 所 
示 。CHIRPS 产品 的 径流 模拟 效果 最 优 (NSE= 
0.72) ,其 次 是 PERSIANN(NSE=0.67) , MSWEP 产品 
的 表现 较 差 (NSE=0.57)。 表 明 CHIRPS 产品 在 黄河 
源 区 有 着 良好 的 适应 性 。 


(a) PERSIANN (b) MSWEP (c) CHIRPS 
~ 259° cc-o98 250 CC=0.95 250r CC=0.98 
om RMSE=9.8 mm T RMSE-14.1 mm T RMSE-8.9 mm 
E 200r kGE=0.94 P 200F kGE=0.95 P KGE=0.95 
i 150 E E 150 E 150 
x Tij Tig 
x z 
m o» k N 
100 xb 100 aj 100 
和 E 3 à 
E: 50} A = 50 E 50 
m = 5 à 
0 50 100 150 200 250 0 50 100 150 200 250 0 50 100 150 200 250 


CMA 降水 量 /mm 月 1 CMA 降水 量 /mm 月 1 


CMA 降水 量 /mm 月 


注 :CC 为 相关 系数 ;RMSE 为 均 方 根 误差 ;KGE 为 Kling-Gupta 系数 。 下 同 。 
图 3 CMA 与 3 套 遥 感 降水 产品 的 降水 散 点 图 


Fig. 3 Precipitation scatter maps of CMA and three remote sensing precipitation products 


(a) ABCD 水 文 模型 率 定期 和 验证 期 的 性 能 


100 € KE >< 验证 >0 

80 
E 60 E 
S Tg 
XS 40 xz 
Kd S 


) V Y D 
1985 


as 
1980 


(b) 不 同和 遥感 降水 产品 的 模拟 性 能 bod 
— Qus — Qsim —— PERSIANN —— MSWEP —— CHIRPS 


100 0 

80 
E i 200 £ 
P Š 
a E | 400 a 

2077 | AAJA AAN UA AN IVA, | | 

0 0 WM OM NJ WM WO W/ NS M OW OW OM WM NM OV CN UV RV V Vv MV oW YY M OM BOV WV V D 600 
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注 :NSE 为 纳什 效率 系数 ;0 为 观测 径流 量 ;Q;, 为 模拟 径流 量 。 
图 4 ABCD 水 文 模型 及 不 同 遥 感 降水 产品 的 模拟 性 能 


Fig.4 Simulation performance of ABCD hydrological model and different remote sensing precipitation products 
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3.3 水 文 干旱 指标 评价 分 析 

1983 一 2018 年 观测 降水 和 遥感 降水 产品 的 SRI 
月 平均 趋势 如 表 3 所 示 ,1983 一 2018 年 观测 降水 和 
各 遥感 产品 的 SRI 均 呈现 增加 趋势 ,但 3 套 遥 感 产 
品 的 SRI 增 加 趋势 要 高 于 观测 降水 ,这 是 由 于 3 套 
遥感 产品 的 原始 降水 数据 的 增加 趋势 均 大 于 观测 
降水 的 结果 。 各 时 间 尺 度 观测 降水 和 遥感 降水 产 
品 的 SRI 均 为 显著 增加 趋势 ,表明 近 36a 源 区 河川 
径流 增加 ,水 文 干旱 趋 缓 。 

1983 一 2018 年 产品 SRI 与 观测 降水 的 SRI 的 相 
关系 数 如 图 $a 所 示 ,3 套 遥感 降水 产品 整体 的 相关 
系数 较 高 , MSWEP 产品 在 各 个 时 间 尺 度 的 相关 系 
数 均 为 最 优 。 随 着 时 间 尺 度 增 大 ,PERSIANN 产品 
的 相关 系数 呈 增 大 趋势 ,MSWEP 和 CHIRPS 产品 无 
明显 趋势 。1983 一 2018 年 遥感 降水 产品 与 观测 降 
水 的 SRI 的 均 方 根 误差 如 图 5b 所 示 ,MSWEP 产品 
在 各 个 时 间 尺 度 的 均 方 根 误差 均 最 小 ， iid 
MSWEP 产 品 与 观测 降水 SRI 之 间 的 偏差 小 ,其 次 是 
CHIRPS 产品 ,PERSIANN 产品 表现 最 差 。1983 一 
2018 年 遥感 降水 产品 与 观测 降水 的 SRI ÉY Kling- 
Gupta 系数 如 图 5c 所 示 ,3 套 产品 的 Kling-Gupta A 


(a) CC (b) RMSE 
PERSIANN MSWEP CHIRPS CC 
0.87 


z 8 


SRII2 SRI6 
SRII2 SRI6 


(d) DEN 


(e) MDD 
CMA PERSIANN MSWEP CHIRPS DEN 


20 
名 
17 
z 
14 & 
11 
8 
5 


SRII 


SRI3 


SRI12 SRI6 
SRI12 SRI6 


注 :DEN 为 干旱 事件 数 


PERSIANN MSWEP CHIRPS RMSE 
0.67 


CMA PERSIANN MSWEP CHIRPS MDD 
15.2 


数 整体 效果 较 差 ,但 PERSIANN 产品 在 年 尺度 
(SRI12) 最 优 。 综 上 ,MSWEP 产 品 在 各 个 时 间 尺 度 
具有 高 相关 系数 , 低 均 方 根 误差 的 特点 ,PERSIANN 
产品 在 SRI12 具 有 高 Kling-Gupta 系数 的 表现 。 

在 较 短 时 间 尺 度 (SRI1 .SRI3) 上 3 套 遥 感 降水 
产品 呈现 显著 低估 (图 $d)。 随 着 时 间 尺 度 增长 ,在 
较 长 时 间 尺 度 (SRI6 .SRI12) 上 3 套 遥 感 降 水 产品 与 
观测 降水 识别 的 干旱 事件 数目 较为 接近 ,表明 遥感 
降水 产品 在 长 时 间 尺 度 能 够 较 好 地 捕捉 干旱 事 
件 。 在 较 短 时 间 尺 度 (SRI1 .SRI3) 上 3 套 遥 感 降 水 


表 3 观测 数据 及 遥感 降水 产品 SRI 趋 势 
Tab.3 SRI trend of observation data and remote sensing 


precipitation products 


不 同时 间 遥感 降水 产品 

尺度 SRI idis PERSIANN — MSWEP CHIRPS 
SRII 0.0016 ^ 0.0042" 0.0034 0.0035" 
SRI3 0.0016" 0.0044" 0.0035 0.0036 ^ 
SRI6 0.0016” 0.0043 ^ 0.0035 0.0037 
SRI12 0.0015" 0.0040 * 0.0035 * 0.0036 ^ 


注 :**## 表 示 通 过 99% 显 著 性 检验 ;SRI1 .SRI3 .SRI6 .SRI12 分 别 表 
示 1 个 月 .3 个 月 .6 个 月 12 个 月 4 个 时 间 尺 度 的 SRI。 


(c) KGE 
PERSIANN MSWEP CHIRPS KOR ai 


SRII 


SRI3 


SRII2 SRI6 


—12.45 


(f) MDS 


CMA PERSIANN MSWEP CHIRPS MDS, r 


12.81 18.71 


e 
SRI3  SRII 


9.97 15.22 
8.03 


5.64 


SRII2 SRI6 


3.25 


;MDD 为 平均 干旱 历时 ;MDS 为 平均 干旱 烈度 ;SRI1 SRI3 SRI6 ,SRII2 分 别 


表示 1 个 月 .3 个 月 .6 个 月 .12 个 月 4 个 时 间 尺 度 的 SRI。 


图 5 通 感 降水 产品 捕捉 水 文 干旱 特征 


Fig. 5 Hydrological drought characteristics captured by remote sensing precipitation products 
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T d EL E Th (Al Se). ZEB AY ISIN HE 
(SRI6) EMSWEP 产 品 与 观测 降水 的 平均 干旱 历时 
最 为 接近 ,PERSIANN 和 CHIRPS 产品 呈现 高 估 倾 
向 。 在 年 尺度 (SRI12) E PERSIANN 产品 与 观测 降 
水 的 平均 干旱 历时 最 接近 , MSWEP 和 CHIRPS 产品 
均 呈 现 显著 低 估 。 在 较 短 时 间 尺 度 (SRI1 .SRI3 ) 上 
3 套 遥 感 降水 产品 均 呈 现 显著 高 估 ( 图 $0) 。 在 较 长 
时 间 尺 度 (SRI6) 上 MSWEP 产品 呈现 略微 高 估 ， 
PERSIANN 和 CHIRPS 产品 则 呈现 显著 高 估 。 在 年 
尺度 (SRI12) 上 PERSIANN 产品 与 观测 降水 的 平均 
干旱 烈度 最 接近 ,MSWEP 和 CHIRPS 产品 均 呈 现 显 
著 低 估 。 综 上 ,3 套 遥 感 产 品 在 短 时 间 尺 度 (SRII1、 
SRI3) 上 均 高 估 了 干旱 事件 的 历时 、 烈 度 , MSWEP 
产品 在 SRI1 .SRI3 和 SRI6 的 平均 干旱 历时 和 烈度 
相对 于 其 他 2 套 产 品 效果 较 优 。PERSIANN 产品 在 
年 尺度 (SRI12) 上 的 平均 干旱 历时 和 烈度 有 较 优 的 
表现 。CHIRPS 产品 捕捉 水 文 干旱 事件 特征 的 性 能 
不 如 PERSIANN 和 MSWEP 产 品 。 


4 讨论 


诸多 学 者 在 黄河 源 区 开展 相关 研究 。 韩 兰 英 
SE BESS A 1961 年 以 来 黄河 流域 整体 趋 于 暖 干 化 ， 
气候 变化 使 得 径流 呈 减 少 趋势 ,从 而 使 干旱 风险 加 
剧 。 但 本 研究 发 现 1983 一 2018 年 黄河 源 区 观测 降 
水 和 遥感 降水 产品 计算 的 SRI 在 各 个 时 间 尺 度 均 呈 
现 显著 增加 趋势 ,表明 近 36a 河 川 径流 增加 ,水 文 干 
PAM, 5 WE CEO WEE 1991 一 2018 年 黄河 
源 区 的 径流 量 表现 出 增加 趋势 的 结论 一 致 。 

诸多 学 者 评估 了 遥感 降水 产品 在 黄河 源 区 的 
PERE. Bai? 分析 了 5 套 遥 感 产 品 在 三 江 源 区 的 
水 文 性 能 ,得 出 MSWEP 产 品 在 黄河 源 区 优 于 其 他 降 
水 产品 的 结论 ,这 与 本 人 研究 比较 一 怪 。 张 春 知 等 
发 现 PERSIANN 产品 在 黄河 源 区 极 大 的 高 估 了 降水 
fi ,本 研究 发 现 PERSIANN 产品 计算 的 SRI 增 加 趋 
势 约 为 观测 降水 模拟 SRI 的 3 倍 ,间接 说 明了 PERST- 
ANN 产品 在 黄河 源 区 高 估 了 降水 量 。 许 昕 形 等 3 
研究 表明 ,1981 一 2020 年 MSWEP 产 品 在 黄河 源 区 
的 SPI3 倾向 于 低估 干旱 历时 、 烈 度 和 极 值 ,而 本 研 
Zt A Jil 1983—2018 年 MSWEP 产品 在 黄河 源 区 的 
SRI3 倾向 于 高 估 王 旱 历 时 和 烈度 。 本 人 研究 尝试 将 
研究 时 段 扩 展 到 1981—2020 4E. , KHL MSWEP 产品 
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在 2017 一 2020 年 于 黄河 源 区 有 严重 的 高 估 降 水 ,这 
会 给 干旱 事件 的 分 析 带 来 极 大 的 不 确定 性 。 本 研 
究 中 遥感 降水 产品 对 黄河 源 区 均 高 估 了 降水 ,低估 
了 和 干旱 事件 数目 ,高 估 SRIL 和 SRI3 的 干旱 历时 和 
烈度 ,根据 二 等 的 的 研究 ,复杂 的 地 形 和 稀 玻 的 站 
点 会 制约 遥感 降水 产品 的 精度 ,从 而 造成 误差 。 
虽然 本 研究 对 数据 源 进行 严格 的 质量 控制 ,但 
不 可 避免 的 存在 不 确定 性 。 首 先 , 选 取 的 Anuspline 
插值 方法 虽然 考虑 了 高 程 和 地 形 等 影响 ,但 依然 存 
在 误差 。 其 次 ,只 采用 单一 水 文 模型 和 水 文 干旱 
指数 进行 研究 ,可 能 存在 不 确定 性 ,未 来 开展 的 研 
究 中 会 综合 考虑 多 个 水 文 模型 和 水 文 干旱 指数 进行 
对 比 研究 ,以 便 更 加 客观 地 对 遥感 降水 产品 进行 评 
估 。 最 后 ,本 研究 侧重 于 客观 评估 3 套 降 水 产品 对 
于 水 文 干 旱 的 模拟 特征 分 析 ,并 未 定量 订正 3 套 遥 
感 产 品 ,未 来 会 考虑 遥感 产品 对 降水 呈现 高 佑 的 原 
,采取 贝 叶 斯 卷 积 神经 网 络 中 和 机 器 学 习 算法 太 
对 产品 进行 定量 订正 。 


5 结论 


(1) 3 套 遥 感 降水 产品 与 观测 降水 的 时 空 分 布 
格局 整体 一 致 , 其 中 CHIRPS 产品 表现 最 优 。 基 于 
1983—2018 年 观测 降水 率 定 的 水 文 模型 参数 ,将 3 
套 遥 感 降 水 产品 驱动 ABCD 水文 模型 ,CHIRPS 产品 
的 NSE 在 3 套 产 品 中 表现 最 优 ,表明 CHIRPS 产 品 在 
黄河 源 区 有 着 良好 的 适应 性 。 

(2) 1983 一 2018 年 黄河 源 区 各 时 间 尺 度 观 测 降 
水 和 遥感 降水 产品 的 SRI 均 为 显著 增加 趋势 ,表明 
近 36a 源 区 河川 径流 增加 \ 水 文 干 旱 趋 缓 ,但 3 套 遥 
感 降 水 产品 均 呈 现 高 佑 倾向。3 套 遥感 降水 的 SRI 
统计 性 能 存在 差异 ,其 中 MSWEP 产 品 计算 的 SRI 与 
CMA 的 最 为 一 致 ,具有 高 相关 系数 , 低 均 方 根 误差 
和 低 平 均 绝 对 误差 的 特点 。 虽 然 3 套 遥感 降水 产品 
的 Kling-Gupta 系数 总 体 不 高 ,但 在 年 尺度 (SRI12) 
上 PERSIANN 产品 表现 较 优 。 

(3) 1983 一 2018 年 黄河 源 区 的 干旱 事件 根据 游 
程 理论 识别 。3 套 遥感 降水 产品 均 低 估 了 SRI1 和 
SRI3 的 干旱 事件 数目 ,高 估 了 干旱 历时 和 烈度 。 
MSWEP 产 品 对 SRI6 的 模拟 性 能 最 优 ; PERISANN 
产品 在 年 尺度 (SRI12) 与 观测 降水 的 历时 和 烈度 较 
为 接近 。CHIRPS 产品 捕捉 水 文 干旱 事件 特征 的 性 
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能 不 如 PERSIANN 和 MSWEP 产 品 。 
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Simulation performance of remote sensing precipitation products on 
hydrological drought characteristics in the source region 
of the Yellow River 
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Abstract: In regions with scarce data, remote sensing precipitation products provide crucial data for the develop- 
ment of the hydrometeorological disaster mechanism and early warning studies. However, the performance of var- 
ious remote sensing precipitation products exhibits regional heterogeneity. Comprehensively evaluating the per- 
formance of remote sensing precipitation products is critical for their use in hydrometeorological-related research 
and application. Based on this assumption, the study investigated the source region of the Yellow River of China 
by using the observed precipitation data (CMA) from 1983 to 2018 to drive and calibrate the ABCD hydrological 
model. Furthermore, the standardized runoff index (SRI) was used to evaluate the simulation performance of 
three sets of typical remote sensing precipitation products (PERSIANN-CDR, CHIRPS v2.0, MSWEP v2.0) on 
hydrological drought. Furthermore, hydrological drought events were identified by using run theory, and the po- 
tency of remote sensing precipitation to capture hydrological drought characteristics was investigated. The results 
revealed that: (1) The three precipitation products can accurately capture the temporal and spatial distribution pat- 
tern of CMA' s multiyear mean value. Furthermore, the CHIRPS product (Nash-Sutcliffe efficiency coefficient, 
NSE=0.72) outperformed other two products in term of hydrological simulation. (2) The SRI values (SRII, SRI3, 
SRI6, and SRI12) of the four scales simulated by CMA and three sets of remote sensing precipitation products re- 
vealed a significant increase trend (P«0.01), which indicated that the river runoff in the source region increased in 
the last 36 years, and the hydrological drought slowed down. However, the SRI values of the three sets of remote 
sensing precipitation products were overestimated. This result revealed that the deviation correction of precipita- 
tion products in the source area of the Yellow River is necessary. In terms of basic statistical indicators, the SRI 
calculated by the MSWEP product was the most consistent with CMA and was the best. However, on an annual 
scale (SRI12), the PERSIANN product achieved the best performance. (3) Three sets of remote sensing precipita- 
tion products overestimated the drought duration and intensity of SRI1 and SRI3; the MSWEP product achieved 
the best simulation performance for SRI6; and the PERSIANN product has the best simulation performance for 
SRI12. The results of this study can provide scientific decision support for the selection of precipitation product 
data for studying hydrological drought in the source region of the Yellow River. 

Key words: remote sensing precipitation; hydrological drought; ABCD hydrological model; standardized run- 
off index (SRI); the source region of the Yellow River 


